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ABSTRACT 

Using the OVRO, Nobeyama, and IRAM mm-arrays, we searched for "disk"- 
outflow systems in three high-mass (proto)star forming regions: G 16.59—0.05, 
G 23.01-0.41, and G 28.87+0.07. These were selected from a sample of NHg 
cores (Codella, Testi & Cesaroni) associated with OH and H2O maser emission 
(Foster & Caswell) and with no or very faint continuum emission. Our imaging 
of molecular line (including rotational transitions of CH3CN) and 3 mm dust 
continuum emission revealed that these are compact (^ 0.05 - 0.3 pc), massive (~ 
100 - 400 Mq), and hot (~100 K) molecular cores (HMCs), that is hkely sites of 
high-mass star formation prior to the appearance of ultracompact Hii regions. All 
three sources turn out to be associated with molecular outflows from ^^CO and/or 
HCO"'" J =1-0 line imaging. In addition, velocity gradients of 10 - 100 km s~^ 
pc~^ in the innermost (^0.03 - 0.13 pc), densest regions of the G 23.01—0.41 
and G 28.87+0.07 HMCs are identified along directions roughly perpendicular 
to the axes of the corresponding outfiows. All the results suggest that these 
cores might be rotating about the outflow axis, although the contribution of 
rotation to gravitational equilibrium of the HMCs appears to be negligible. Our 
analysis indicates that the 3 HMCs are close to virial equilibrium due to turbulent 
pressure support. Comparison with other similar objects where rotating toroids 
have been identifled so far shows that in our case rotation appears to be much less 
prominent; this can be explained by the combined effect of unfavorable projection, 
large distance, and limited angular resolution with the current interferometers. 

Subject headings: ISM: evolution — ISM: individual (G 16.59-0.05, G 23.01-0.41, 
G 28.87+0.07) — stars: early type — radio continuum: ISM 
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1. Introduction 

The role of disks in the formation process of low-mass stars (i.e. stars with masses <^ 1 
Mq) has been extensively studied in the last two decades through high angular resolution 
observations at various wavelengths. Images of such disks have been obtained in the optical 
(e.g. Burrows et al. 1996) and at mm- wavelengths (e.g. Simon et al. 2000). The latter have 
demonstrated that the majority of the disks undergo Keplerian rotation. These findings are 
consistent with the fact that low-mass stars form through accretion while (partial) conserva- 
tion of angular momentum during the dynamical collapse produces a flattened and rotating 
structure at the center of the core. 

What about high-mass (M* ^ 8 Mq) stars? In this case, formation through accretion 
faces the problem that stars more massive than ~8 Mq reach the zero age main sequence 
still deeply embedded in their parental cores (Palla & Stabler 1993). At this point radiation 
pressure from the newly formed early-type star can halt the infall, thus preventing further 
growth of the stellar mass. Various solutions have been proposed to solve this problem: 
(i) massive stars might form through merging of lower mass stars (Bonnell, Bate & Zinnecker, 
H. 1998; Bonnell & Bate 2002; Bally & Zinnecker 2005); (ii) sufficiently large accretion rates 
could allow the ram pressure of the infalling material to overcome the radiation pressure 
form the star (Behrend & Maeder 2001; McKee & Tan 2003; Bonnell, Vine, & Bate 2004) 
(iii) non-spherical accretion could weaken the effect of radiation pressure by allowing part 
of the photons to escape through evacuated regions along the outflow axis and, at the same 
time, enhancing the ram pressure of the accreting material by focusing it through the disk 
plane (Yorke & Sonnhalter 2002; Krumholz, McKee, & Klein 2005). 

An important test to discriminate between the different hypotheses is the presence of 
rotating, circumstellar disks, which would lend support to the third scenario depicted above. 
This can be determined by inspecting the velocity field of the innermost parts of molecular 
cores where young massive (proto) stars are believed to form. One possibility is to look for 
velocity gradients perpendicular to the direction of the larger scale outflows associated with 
such cores; this would suggest that one is observing rotation about the outflow axis. This 
technique has been adopted by us and other authors successfully inferring the existence of 
rotation in the gas enshrouding high-mass young stellar objects (YSOs) and leading to the 
discovery of circumstellar disk-like structure or "toroids" in a limited number of cases (see 
Cesaroni et al. 2007 for a review on this topic) . 

The goal of the present study was to establish whether the presence of rotation is 
common in high-mass star forming regions by observing three more objects, expected to 
be sites of deeply embedded OB (proto)stars. The ideal target for this type of studies are 
hot molecular cores (HMCs), which are believed to be the cradles of massive stars (see e.g. 
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Kurtz et al. 2000; Cesaroni et al. 2007). Beside other studies, the one by Codella, Testi 
& Cesaroni (1997, hereafter CTC97), has proved successful in identifying the birthplaces of 
young massive stars as dense NH3 cores associated with OH and H2O maser emission. One 
of these HMCs (G 24.78+0.08) has been the subject of a series of articles by us (Furuya et 
al. 2002; Cesaroni et al. 2003; Beltran et al. 2004, 2005), which have shown that this is 
a unique object, characterized by the simultaneous presence of rotation, outflow, and infall 
towards a hypercompact Hii region ionized by an 09.5 star (Beltran et al. 2006). With this 
in mind, we completed the study of the sources in CTC97 by observing three more objects 
in the same tracers used to investigate G 24.78+0.08. The selected targets for this study are 
G 16.59-0.05 atd = 4.7 kpc, G 28.87+0.07 at 7.4 kpc, and G 23.01-0.41 at 10.7 kpc. The 
fact that in all three cases no or only faint free-free continuum emission has been detected, 
notwithstanding the large luminosities of the sources (CTC97), suggests that the embedded 
YSOs might be in an even earlier evolutionary phase than those in G 24.78+0.08 (Furuya et 
al. 2002; Beltran et al. 2004). 



2. Observations and Data Retrieval 

2.1. Millimeter Array Observations 

Aperture synthesis observations of molecular lines and continuum emission at 3 mm 
were carried out using the Owens Valley Radio Observatory (OVRO)E| millimeter array to- 
wards G 16.59-0.05 and G 28.87+0.07, and the Nobeyama Millimeter Array (NMA) of the 
Nobeyama Radio Observatorj]^ towards G 23.01—0.41. We also carried out ^^CO (1-0) line 
and millimeter continuum emission imaging with the Plateau de Bure Interferometer (PdBI) 
of the Institut de Radio Astronomie Millimetriqu^ (IRAM). The parameters for the contin- 
uum and molecular line observations are summarized in Tables [H and [21 respectively. 



^Research at the Owens Valley Radio Observatory is supported by the National Science Foundation 
through NSF grant AST 02-28955. 

^Nobeyama Radio Observatory is a branch of the National Astronomical Observatory, operated by the 
Ministry of Education, Culture, Sports, Science and Technology, Japan. 

3IRAM is supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain). 
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2.1.1. OVRO Observations 

The OVRO observations of G 16.59— 0.05 and G 28.87+0.07 were carried out in the 
period from 2003 September to 2004 May in 3 array configurations (E, H, and UH). The 
shortest projected basehne length, i.e. the shadowing limit, is about 12.8 m. This makes our 
OVRO observations insensitive to structures more extended than 54", corresponding to 1.2 
and 1.9 pc at the distances of G 16.59— 0.05 and G 28.87+0.07, respectively. We observed 
the GH3GN (5-4) and N2H+ (1-0) hues in the upper sideband (USB) and the HGO+ (1-0) 
line in the lower sideband (LSB). For the continuum emission, we simultaneously used the 
Gontinuum Gorrelator with an effective bandwidth of 4 GHz and the newly installed GOBRA 
with 8 GHz bandwidth. We configured the digital correlator with 31 MHz bandwidth and 62 
channels centered at 91979.970 MHz, so to cover the X = to 3 components of the GH3GN 
(5-4) transition. We used 3G273 and 3G 454.3 as passband calibrators, and NRAO530 and 
J1743-038 as phase and gain calibrators respectively for G 16.59-0.05 and G 28.87+0.07. 
The flux densities of NRAO 530 and J1743— 038 were determined from observations of Uranus 
and Neptune. We estimate the uncertainty of flux cahbrations to be 10%. The data were 
cahbrated and edited using the MM A and MIRIAD packages. We constructed continuum 
images from the GOBRA data by employing the multi-frequency synthesis method: the flnal 
center frequency is 90.584 GHz. 



2.1.2. NMA Observations 

Our NMA observations towards G 23.01—0.41 were carried out in the period from 2003 
December to 2004 May with 3 array configurations (D, G, and AB). The largest detectable 
source size is 50", corresponding to 2.6 pc at c? = 10.7 kpc. We observed the GH3GN (6-5), 
HNCO (5-4), ^^CO (1-0) and C^^O (1-0) lines in the USB. For the continuum and fine 
emission, except for CH3CN, we employed the Ultra Wide Band Correlator (UWBC) with 
a 512 MHz bandwidth in each sideband; this configuration gives a total bandwidth of 1 GHz 
for the continuum and a velocity resolution of 5.5 km s^^ for the lines. We used the FX 
correlator with 32 MHz bandwidth centered at 110374.0 MHz, thus covering the X = to 3 
components of the GH3GN (6-5) transition. We adopted 3G 273 as passband cahbrator and 
J1743— 038 as phase and gain calibrator. The fiux densities of J1743— 038 were bootstrapped 
from Uranus, and the uncertainty of the flux calibration is estimated to be 10%. We decided 
not to use the ^^CO and G^®0 data taken with the G and AB configurations as these are 
too extended and resolve the emission. All the data were calibrated and edited using the 
UVPR0G2 and MlRlAD packages. 
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2.1.3. PdBI Observations 

Our ^^CO (1-0) and 2.6 mm continuum emission observations towards the 3 HMCs were 
carried out with the IRAM 5-antenna interferometer on Plateau de Bure in 1998 April and 
May. The D and CI configurations were used, yielding a largest detectable angular scale 
of ~44". The 82 - 116 GHz SIS receivers were tuned in SSB at the frequency of the ^^CO 
(1-0) line, and the facility correlator was configured with a bandwidth of 40 MHz centered 
at the same frequency. Continuum emission has also been measured with two bandwidths of 
160 MHz each. Line free channels were averaged to produce a continuum image which was 
then subtracted from the line data in the visibility plane. The fiux densities of J1833— 210 
(1.7 Jy) and J1743— 038 (3.0 Jy) were bootstrapped from 3C273 whose fiux density was 
assumed to be 16.7 Jy. We estimated an overall uncertainty on fiux calibration of 20%. 
Data editing, calibrations, and image construction were done using the GILDAS software 
package developed at IRAM. Tables [1] and [2] summarize the observing parameters for the 
PdBI continuum and line imaging, respectively. 

2.2. GLIMPSE and MSX Archive Data 

We have retrieved infrared (IR) images at 21 fim from the Midcourse Space Experiment 
(MSX; Price et al. 2001) survey and at 8.0, 5.8, 4.5, and 3.6 /xm from the Galactic Legacy 
Infrared Mid-Plane Survey (GLIMPSE) survey making use of the Infrared Array Camera 
(IRAC) on-board the Spitzer satellite. The calibrated data from the Spitzer Science Center 
were processed through the GLIMPSE pipeline reduction system (Benjamin et al. 2003; 
Whitney et al. 2004). All the infrared (IR) images were used without spatial smoothing; 
the pixel size is 6" and ~ 1'.'2 for the MSX and IRAC data, respectively. 

3. Results and Discussion 
3.1. Continuum Emission 

3.1.1. Source Identification Based on Millimeter Images 

Figure [1] compares our 3 mm continuum emission maps to the IR images. The latter 
correspond to the MSX data at 21 /xm and the GLIMPSE ones at the longest (8.0 /im) and 
shortest (3.6 yum) wavelengths. Here we do not show GLIMPSE images at 5.8 and 4.5 /xm 
as they are basically similar to those at 8.0 and 3.6 ^m. All 3 objects show intense, compact 
3 mm continuum emission coincident with IR emission, although one must consider that the 
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angular resolution of the MSX images is much worse than that of our interferometric maps. 
In the miUimeter maps of G 28.87+0.07, one can also see weaker 3.3 mm continuum sources 
located ~ 14" to the north (hereafter G28.87+0.07B and C; see Figure [T}d) and a 2.6 mm 
source ~ 10" to the east (G 28.87+0.07D; same figure) of the main core found by CTC97. 
No other mm continuum emission is detected with S/N > 5 in our OVRO, NMA, and PdBI 
fields of view towards the 3 sources. Assuming that the flux density Su-, varies as z/(^+^) with 
(3 = 1.5 (Preibisch 1993; Molinari et al. 2000), the expected flux of G 28.87+0.07 B and C 
in the PdBI beam at 2.6 mm is below the sensitivity of our images (3.5 mJy corresponding 
to 2.2o"). Similarly, assuming that G 28.87+0.07 D matches the PdBI beam, this source is 
resolved by the OVRO beam and its flux at 3.3 mm corresponds to only lAa of the OVRO 
continuum image, thus explaining why it is not detected. 

The GLIMPSE images of G 16.59-0.05 and G 28.87+0.07 show IR emission both from 
the main cores and from their surroundings. On the other hand, G 23.01—0.41 seems rather 
isolated, although we may be missing weaker objects due to the large distance to the region. 
In the following we make no attempt to identify other IR sources beside the counterparts 
of those detected at mm wavelengths, because this would require an analysis which goes 
beyond the purposes of the present study. Notice that no other mm sources with 3< S/N < 
5 is seen in the IRAC bands. We summarized the identified continuum sources in Table [3l 
where peak positions and effective radius (-Reff) of the mm-continuum sources as well as their 
S,y integrated inside the 5a level contours are given. Here the observed angular radius of the 
source is computed from \/ Ajix where A is the area enclosed by the 50% contour level of 
the emission. Deconvolution assuming the source and the synthesized beam to be Gaussian 
is then applied, thus obtaining i?eff) assuming that the emission has a Gaussian distribution. 
It should be noted that i?eff does not depend on the image noise levels. To estimate the core 
mass from the continuum flux the latter has been integrated all over the emitting region, 
namely inside the corresponding 5cr contour level. The fluxes of the IR counterparts are 
listed in Table |H 

3.1.2. Continuum Spectra and Core Masses Estimated from the 3 mm Flux Densities 

Figure [2] shows the continuum spectra of the 3 HMCs, as from Tables [3] and |H No 
contamination by free-free emission is expected to affect our 3 mm fluxes, as the upper 
limits obtained by CTC97 at 1.26 cm guarantee a maximum free-free flux of ~ 1.4 mJy at 
3.0 mm - assuming optically thin emission. Establishing the origin of the IR emission requires 
a model fitting of the spectral energy distribution (SED) that goes beyond the purpose of 
this study. We only note that a simple grey-body fit cannot consistently reproduce the SED 
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from the mm to the near-IR regime. This suggests that the structure of the cores is rather 
complex, probably hosting a multiple stellar system such as the one seen in the G 29.96—0.02 
star forming region: here, recent sub-arcsecond resolution (~ 1800 AU resolution in linear 
scale) Submillimeter Array (SMA) imaging resolved the HMC (Maxia et al. 2001; Olmi et 
al. 2003) into six submm continuum sources (Beuther et al. 2007). Although our angular 
resolutions in AU are 5-12 times worse, by analogy to the case of G 29.96—0.02, we may 
argue that multiple sources are present in our cores, which might explain why the SEDs 
cannot be fitted with a simple grey-body fit. Furthermore, temperature gradients as well as 
optical depth effects and dumpiness should be taken into account to obtain a satisfactory 
fit. In particular images with arcsecond resolution at wavelengths above 20 /im would be of 
great importance to obtain a precise estimate of the bolometric luminosity (and hence of the 
mass) of the embedded (proto) stars. 

The 3 mm flux density in Table [3] can be used to estimate the mass of the cores (Mdust), 
under a few assumptions on the dust emissivity. More specifically, we use Hildebrand's 
(1983) relation Mdust= ^.g^rLt) ^^^^ = ^'^ ^ 10-='(z//231 GHz)^ cm^ g-i and (3 = 1.5. 
When calculating Mdust, we set Tdust equal to the rotational temperature obtained from 
the CH3CN lines (see §3.5.11) . assuming that gas and dust are well-coupled (e.g., Kriigel & 
Walmsley 1984), a reasonable hypothesis at densities as high as those traced by the CH3CN 
emission (typically > 10^"'' cm~^). 

The derived values of Mdust range from 95 to 380 Mq (Tabled]). Our 3 targets appear 
to fall in the mass range for "heavy" HMCs 100 Mq] Cesaroni et al. 2007). Cores that 
massive are likely to host multiple stars, as opposed to "light" HMCs (< a few 10 Mq), 
which appear to contain only single OB stars or binary systems (Cesaroni et al. 2007). 

In G 28.87-1-0.07 we have identified 3 additional weaker sources. Assuming that Tdust 
is equal to the excitation temperature of N2H"'"(l-0) (see §3.3.21 and Figure Et), we obtain 
Mdust= 7M0 with Tdust= 46 K for G 28. 87-1-0. 07D. The excitation temperature is estimated 
from the hyperfine structure analysis of mean spectra of the N2H"'" emission towards a con- 
densation hosting G28.87+0.07D. If we assume that Tdust^ 46K is valid for G28.87+0.07B 
and C, we obtain Mdust of 29 Mq for G28.87+0.07B and 28 Mq for G 28.87+0.07C. These 
results suggest that the G 28.87+0.07 HMC (Mdust= lOOM©; Table E]) is surrounded by less 
massive objects. Such a situation appears to be typical of high-mass star forming regions, 
where embedded OB (proto) stars happen to be surrounded by less embedded, possibly more 
evolved, lower mass stars (Molinari et al. 1998; Fontani et al. 2004a, 2004b). Further- 
more, because of the absence of free-free emission (CTC97), we argue that the G 28.87-1-0.07 
cluster is much younger than other clusters of massive YSOs, such as, e.g., the cluster of 
ultracompact Hii regions in G 19.61—0.23 (Furuya et al. 2005). 
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Finally, in order to estimate the CH3CN abundance in the cores (see §3.5.2p . we have 
calculated the mean molecular hydrogen column density, (A^(H2)), over the cores. To make 
the comparison between continuum and CH3CN emission as consistent as possible, we re- 
constructed the 3 mm images by fixing the beam sizes to those of the CH3CN images (Table 
[2]). We subsequently integrated the continuum fiux over the region enclosed by the 5a level 
contours. In Table [5] we summarize the calculated mean column densities ((A^(H2))) which 
are (3 - 4)xl0^^ cm~^, implying optical depths on the order of 10"'^, consistent with our 
assumption that the dust emission is optically thin at 3 mm. 



3.2. Spectral Line Profiles towards the Cores 

Figure [3] shows the line spectra at the peak positions of each integrated intensity map, 
in units of brightness temperature in the synthesized beam (Tgb). For the CH3CN emission, 
we show the K = 3 component as this is not blended with the other iiT-ladder lines. 

G 16.59—0.06 The ^^CO and HCO"*" peak spectra present deep "absorption" features 
around the systemic velocity, the latter being estimated from the CH3CN lines (see §3.5.ip . 
Most likely this is not due to real self-absorption but to the fact that emission at low veloc- 
ities arises from extended regions resolved out by the interferometer, as found in the case of 
G24.78+0.08 (see Cesaroni et al. 2003). However, only single-dish mapping of the regions 
under study in the same lines can establish the nature of such absorption. 

The ^^CO spectrum shows prominent high- velocity wings on both sides of the line, 
indicating the presence of a molecular outfiow. In order to adopt suitable velocity ranges for 
the outfiow lobes traced by the ^^CO emission, we define terminal (Vt) and boundary (Vb) 
velocities: Vt is the most blue- and redshifted LSR-velocity where the wing emission drops 
below the 1.5a level and Vb's are the velocities where the intensity of the CH3CN K = 3 line 
falls below 1.5a (see caption of Figure H] for the Vt and Vh values). 

G 28.87+0.08 The ^^CO and HCO+ (1-0) spectra towards G 28.87+0.07 are very similar 
each other in the sense that the emission around V^ys is almost missing, and the blueshifted 
wing is barely detected, whereas prominent redshifted emission is seen. The CH3CN K = 3 
emission peaks around Vlsr — 103 km s~^ where the ^^CO and HCO"'" spectra present a 
sharp edge. Also in this case, as for G 16.59—0.05, we argue that the spectral profiles of the 
^^CO and HCO+ lines may suffer of missing short spacing in our interferometric observations. 
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G 23.01-0.41 The profile of the ^^CO fine towards G 23.01-0.41 presents a redshifted 
wing much more prominent than the blue one, while the emission close to Vsys is strongly 
depressed. On the other hand, the lines of the lower density tracers ^'^CO and C^^O present 
well defined peaks, although the velocity resolutions are poor. The peak LSR-velocities of the 
^^CO and C^^O emission shifted slightly towards the blue with respect to the peak velocity 
of CH3CN - which we have assumed as representative of the bulk material. 

3.3. Maps of the Molecular Line Emission 

3.3.1. G 16.59-0.06 

Figure H] compares maps of molecular line and continuum emission in G 16.59—0.05. 
Panel (a) shows an overlay of the ^^CO wing emission on the 3 mm continuum emission. The 
wing emission is integrated over the velocity ranges between Vb and Vt for both the blue 
and red lobes ( §3.21) . The G16.59— 0.05 outfiow map shows a clear bipolarity along a line 
with P.A.~135°. We define the outfiow axis as the line connecting the peaks of the blue lobe 
and 3 mm continuum map. The blue lobe shows a remarkable structure with a large opening 
angle of ~ 180°, whereas the red lobe is much more coUimated. Interestingly, both lobes have 
a large aspect ratio (> 1), defined as the ratio between the width and the length of the lobe. 
Moreover, the two lobes have a rather large spatial overlap, coinciding also with the position 
of the 3 mm continuum emission. Such a geometry indicates that the inclination angle of the 
outfiow axis with respect to the line-of-sight - albeit difficult to quantify precisely - must be 
rather small. 

In Figure Hi, we see that, while the distributions of CH3CN (5-4) = + 1 and NH3 
(3,3) emission agree well with each other, the peak position of the 3.3 mm source is displaced 
by ^ 1'.'4 to the northwest of the CH3CN peak. We also verified that the peak positions 
of the K = 2 and 3 emission agree with that of the K = + 1 emission. Note that the 
3.3 mm continuum and CH3CN visibilities shared the same calibrators as they were observed 
simultaneously ( §2.1.ip . and that the peak positions of the 3.3 and 2.6 mm continuum sources 
(see Table [3]), measured respectively with OVRO and PdBI, agree with each other within the 
errors. These facts indicate that the positional displacement between the CH3CN and 3 mm 
continuum emission is real (see §3.61 for a discussion) . One cannot rule out the possibility 
that such a displacement is due to optical depth effects, although this seems unlikely as 
line emission peaks at the same position for all species (NH3 and CII3CN), and the CH3CN 
opacity quite low (^ 2; described in §3.5. ip . Observation of higher-excitation CH3CN lines 
might help shedding light on this issue. The peak position of the more extended HCO^ 
emission agrees with those of NH3 and CH3CN, suggesting that the HCO^ (1^0) line also 
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traces the HMC On the other hand, the N2H"'" emission has a totally different morphology 
from the other molecules, in the sense that it has a fairly extended structure and its peak 
position differs from both those of the other lines and that of the 3 mm continuum. The N2H"'" 
emission is more extended than the HCO"'" emission, suggesting that it might be associated 
not only with the HMC, but also with the other sources seen in the IR images (Figure [1^). 

To estimate the column density of the region, we analyzed the HCO^ and N2H^ data 
in the following fashion. The mean column density of a linear rigid rotor is calculated from 
Eq. (Al) of Scoville et al. (1986). Since we do not have measurements of the optical depth 
(r) and excitation temperature (Tex) of the HCO^ line, we assumed that the HCO"'" emission 
is optically thin and estimated the column density for two extreme values T^y^: a minimum 
Tex, obtained from the line peak Tgb (32.8 K from Figure [3^) and a maximum Tgx, from 
Trot of the CH3CN emission (130 K; described in §3.5.ip . Correspondingly, the mean HCO^ 
column density over a region enclosed by its 5a contour level is A^hco+ — 6.3 x 10^^ cm~^ 
and 2.1 x 10^^ cm~^. For the fractional abundance, X{B.CO~^), we assume a value of 10^^ 
(see e.g. van Dishoeck et al. 1993). Under this assumption, the mass traced by the HCO"'" 
(1-0) line (Mlte) ranges from 440 to 1500 Mq (Table [6]). We stress that these estimates 
are to be taken as lower limits, as part of the line emission close to the systemic velocity is 
resolved out in our interferometric maps. 

Subsequently we analyzed the N2H"'" data because the hyperfine structure (HFS) of 
the transition allows an estimate of the line optical depth (see, e.g., Benson, Caselli & 
Myers 1998; Caselli et al. 2002). The details of our N2H"'" HFS analysis are described 
in Appendix B of Furuya et al. (2006). For this purpose, we used two spectra obtained 
integrating the emission over the 5cr contour level of the 3 mm continuum emission and the 
integrated N2H+ line emission. Note that the N2H"'" emission does not peak at the position of 
the 3 mm continuum source. For the first spectrum, our analysis gives Tex — 41 ± 9 K, total 
optical depth (corresponding to the sum of the optical depths of the 7 hyperfine components) 
Ttot — 0.28 ± 0.03, and intrinsic line width Avi^t — 2.10 ± 0.04 km s~^ towards the 3 mm 
source. From the second spectrum, we obtain Tgx = 16±2 K, Ttot = 0.91 and At>int = 
2.82±0.02 km s"^. Using Eq. (B7) of Furuya et al. (2006), we calculate a column density 
iVN2H+ — (7.8 ± 1.5) X 10^'^ cm~^ towards the 3 mm continuum source and (6.1 ± 0.4) x 10^^ 
cm~^ for the average spectrum. If we assume X(N2H"'") = 3 x 10""^° suitable for high-mass 
star forming regions (Womack et al. 1992), the mass estimate for the entire region is 540 
Mq (Table E]). 

These numbers demonstrate that the N2H"'" line is coming from a lower (column) density 
region, distinct from the HMC traced by the continuum and CH3CN line emission - whose 
analysis we defer to §3.5[ This is especially evident from the line width, which for the N2H"'" 
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line is ~50% of that of the CH3CN transitions (see Table [8]). 

3.3.2. G 28.87+0.06 

The ^^CO spectrum of G 28.87+0.07 strongly suggests the presence of high-velocity 
outflowing gas ( §3.21) . However, the maps of the wing emission presented in Figure Eti' do 
not seem to outline a clear bipolar structure with respect to the 3 mm source position. A 
prominent red lobe lies to the southwest, whereas the blue lobe shows an elongated structure 
aligned approximately east- west as well as a condensation to the northeast of the main core. 
If we assume that all the CO lobe is associated with an outflow powered by a source in the 
core, the observed morphology and velocity structure can be interpreted in two ways: (i) one 
outflow, with patchy emission from the lobes and with a large opening angle (close to 90°) 
around the axis denoted by the solid line in Figure [S^i; (ii) two outflows, well collimated along 
the axes represented by the dashed hues in the same figure. On the basis of the present data 
it is impossible to discriminate between these two hypotheses. Given the existence of several 
sources in G 28.87+0.07 ( §3.1.ip and the fact that multiple outflows appear to be frequent in 
high-mass star forming regions (see Beuther et al. 2002b, 2003, 2004), one cannot rule out 
the possibility that the blueshifted condensation to the northeast and some portion of the 
eastern blueshifted lobe are associated respectively with G28.87+0.07B (and/or C) and D. 
However, for the sake of simplicity, in the following we will consider only the two scenarios 
proposed above, which assume that it is the main G 28.87+0.07 core to host the source(s) 
powering the outflow(s). 

In G 28.87+0.07, the peak positions of all tracers but N2H+ roughly coincide with the 
3 mm continuum peak (FigureOi). Both the NH3 and CH3CN emission are compact, whereas 
the HCO+ emission is slightly elongated in the north-south direction (Figure [5]d). We believe 
that the HCO"^ emission mostly arises from the HMC and not from the outflow lobes because 
its elongation is not parallel to any of the possible outflow axes and only very faint HCO"^ 
line wings are detected (see Figure [3)d). As done for G 16.59—0.05, we estimate iVHCo+ for 
the region enclosed by the 5(T level contour assuming two extreme values for the excitation 
temperature. We obtain A^hco+ = 3.9 x 10^^ cm"^ for Tex = 96 K, and 9.2 x 10^^ cm~^ 
for Tex = 17 K, yielding a mass range Mlte— 84 - 360 Mq (Table [6]). Here, the maximum 
Tex is assumed equal to the value of Trot obtained from the CH3CN emission. On the other 
hand, the minimum Tex is taken equal to the peak T^h of the HCO+ spectrum [14.5 K from 
Figure El Tex- TPf''''(HCO+) + Tbg = 17 K, where Tbg is the cosmic background tempera- 
ture]. Independently of the value chosen, this estimate represents a lower limit because of 
the different UV coverages, beam dilutions, and optical depth effects. We stress that the 
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estimated A^hco+ varies only by a factor of 4.3 regardless of the large uncertainty on Tgx- 

It is interesting that the N2H+ emission in G 28.87+0.07 cannot be seen towards the 
center of the main core, where instead NH3 emission is detected (CTC97). The distribution 
of the N2H"'" emission is very patchy and several condensations are seen close to the edge 
of the ^^CO outflow lobes (Figure (St). The eastern N2H+ condensation seems to be host 
G 28.87+0.07D, although the N2H"'" peak position does not agree with the mm peaks. In 
fact, none of the N2H"'" condensations matches the 3 mm source positions, which may indicate 
a low N2H"'" abundances in these cores, as suggested by Womack et al. (1992), who argue 
that N2H"'" is likely to be destroyed in regions of very high density (riHa > 10^ cm~^) and 
high temperature (Tk> 200 K). Alternatively, one might speculate that N2H+ could trace 
the borders of a region evacuated by the outflow. However, our fit to the N2H+ lines give 
LSR velocities of around 103 km s~^ for these condensations, namely too similar to the V^ys 
to be typical of gas participating in the outfiow motion. 



3.3.3. G23.OI-O.4I 

We have seen that the red wing of the ^^CO spectrum towards G 23.01—0.41 is more 
prominent than the blue wing (Figure [3t). Figure presents an overlay of the high velocity 
^^CO emission map on the 3 mm continuum image; FigureslHt and d, respectively, show such 
overlay maps for the ^'^CO and HNCO wing emissions. To obtain these maps, we integrated 
high velocity emission between the Vb and Vt for each wing (see §3.21 for their definitions). 
In order to keep a consistency in selecting velocity ranges for the 3 lines as much as possible, 
we integrated the blueshifted ^^CO emission emission up to Vlsr = 49.9 km s~^ beyond the 
Vt(blue) of 64.6 km s~^ which is the first LSR-velocity where the blueshifted wing drops below 
the 1.5(7 level (see Figure (St). The resultant velocity ranges for the blueshifted gas are almost 
the same for the 3 tracers, while velocity interval for the ^^CO red wing is approximately 
twice of those for the redshifted ^^CO and HNCO emission. 

In Figure Et, a clear bipolar structure is seen with the bright redshifted lobe lying 
to the southwest and the blueshifted lobes to the northeast and southwest. In addition, 
a faint redshifted emission is found to the northeast of the 3 mm continuum peak. The 
velocity structure traced by the ^^CO and HNCO seems to reconcile with that by the ^^CO 
line in the sense that the bright blueshifted gas lies to the southwest of the 3 mm source, 
although the HNCO does not show blueshifted gas to the northwest. We also point out that 
a fainter redshifted HNCO condensation is seen to the northeast, whose position matches 
the northeastern ^^CO redshifted lobe. In addition, one may notice that the redshifted gas 
outlined by the ^^CO are basically same as the ^^CO; both the emission show their peaks to 
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the south of the 3 mm source and elongated towards southwest and northeast. 

As in the case of G 28.87+0.07 (see §3.3.21) . the above picture of the high velocity 
wings may be interpreted either with the existence of two distinct outflows or with the fact 
that the outflow axis is sufficiently close to the plane of the sky to allow the observation 
of both blue- and redshifted gas in each lobe. In the former hypothesis, one of the two 
putative (proto) stars drives the northeastern blueshifted lobe as well as the southwestern 
redshifted one, and the other (proto)star does the remaining two lobes. If this is the case, 
each (proto)star is driving a pair of outflow lobes whose masses (hence, momentum rates, 
see Table [7]) differs significantly, which seems not to reasonable. In conclusion, we prefer the 
single outfiow hypothesis with the following two supporting reasons: (1) the high velocity 
gas maps outlined by the ^^CO and HNCO do not contradict with the velocity structure 
traced by the ^^CO; neither ^'^CO nor HNCO maps supports the presence of two outfiows, 
and (2) we identified only a 3 mm continuum source towards the center of the HMC ( §3.1.ip . 
We, thus, adopt the solid line in Figure [6^ as the outfiow axis. This line passes across the 
3 mm peak position and as close as possible to the peak positions of the blue- and redshifted 
HNCO lobes. 

In Figures [6)d and c we show maps of the ^^CO and C^^O bulk emission integrated over 
a velocity range shown by the green bars in Figure [3t;. These arise from a region ten times 
larger than the HMC traced by higher density tracers such as CH3CN. In addition, the 
emission appears to peak a few arcsecond to the south of the HMC, which may be explained 
if the ^'^CO and C^^O lines are partially thick, thus hiding the densest part of the molecular 
cloud. To confirm this hypothesis, we have estimated the optical depth of the ^^CO and 
C'^^O gas from the ratio of the two isotopomers, assuming a relative abundance of 5.5. Note 
that our estimate is not affected by uncertainties due to relative calibration errors or different 
UV sampling, because the two lines were observed simultaneously ( §2.1.ip . Although the 
opacity estimate turns out to vary slightly across the line, depending on the velocity (from 
< 1 to 3.4), we conclude that a mean value of r[^^CO (1 — 0)] = 2 is appropriate for our 
purposes. This indicates that the optical depth of the ^'^CO (1-0) line may be sufficiently 
high to prevent the detection of the HMC. On the other hand, C^^O should be less affected 
by this problem. In fact a tail of C^^O (1-0) bulk emission is seen in Figure towards 
the position of the HMC, although the main peak is still shifted to the south. This may be 
also the effect of low angular resolution, which privileges regions with larger beam averaged 
column density with respect to those (like the HMC) with enhanced volume and source 
averaged column density. 

We thus obtain mean column densities of 4.2x10^^ cm~^ for ^'^CO and 7.1x10^^ cm^^ 
for C^^O, over the regions enclosed by the corresponding 5a contour level, which imply Mlte 
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of 1x10^ Mq for ^^CO and 2x10^ Mq for C^^O (see Table [6]). Here we assume that the 
same T^x for both ^'^CO and C^'^O and take it equal to the peak Tsb of the ^^CO emission 
[ 10.3 K from Figure El T^^^ T^^^^{^^CO) + Tbg = 13 K]. Despite the lack of zero-spacing 
information in these interferometric maps, masses like these compare well to those measured 
in similar objects with single-dish telescopes (e.g., G 24.78+0.08; Cesaroni et al. 2003). 

Lastly, in Figure [6li we show the map of the CH3CN emission integrated over the blue 
and red bars in Figure One may see that the shape of the HMC traced by this and 
the NH3 line (Figure Eb), as well as the 3 mm continuum (Figure EtL) is slightly elongated 
perpendicularly to the outflow axis. It is thus tempting to speculate that one is observing a 
disk-like structure rotating about that axis. We will discuss this issue in §3.61 

3.4. Physical Properties of the Molecular Outflows 

In order to calculate the outflow properties, we used the ^^CO wing emission maps shown 
in Figures HHHl For G 28.87-1-0.07, where multiple outflows might be present (see §3.3.2p . we 
consider only the most prominent lobes, as these represent the dominant contribution to the 
mechanical luminosity and momentum rate of the outflow system. A caveat to this approach 
is that the secondary outflow might be weak because heavily resolved in our interferometric 
maps. In other words, one cannot rule out the possibility that better UV sampling at short 
spacing could reveal a stronger, better defined outflow pattern. Only single-dish mapping of 
these regions may shed light on this issue. 

The outflow parameters are calculated as follows. We integrated the outflow lobe emis- 
sion encompassing the 5a level contours to calculate the mass of molecular hydrogen consti- 
tuting the lobes (Miobe), assuming that the wing emission is optically thin. Indeed, we have 
estimated the optical depth in the ^^CO line wings for G 23.01—0.41, for which ^^CO data 
are available. We first reconstructed the NMA ^^CO data with the same CLEAN beam as 
the PdBI ^^CO data and then smoothed the latter to the same spectral resolution as the 
former. From the ratio between the ^^CO and ^^CO line profiles we computed the optical 
depth for the blue (47.2 < VLSR/km s-^< 69.0) and red (79.9 < Visn/km s-^< 96.2) wings, 
which are respectively 0.70 ± 0.2 and 1.1 ± 0.5. This indicates that the correction due to 
optical depth effects for the outflow mass and the relevant parameters is very small, a factor 
~1.5. Given the similarity in the source properties, we assumed that the wing emission is 
optically thin also in G 16.59— 0.05 and G 28.87+0.07 and used the same approach as for 
G23.01-0.41. 

To calculate the outflow parameters, as a lower limit for Tgx we adopted the maximum 
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Tsb observed in all low-density tracers (CO and HCO"'"; see Figure[3]), namely Tex = 35K for 
G 16.59-0.05, 31 K for G 28.87+0.07, and 15 K for G 23.01-0.41. Subsequently we calculated 
the kinematical properties of the flow (see Table [7j) such as the dynamical time scale (td), 
mass loss rate (Mfjow); and momentum rate (-Fco)- The uncertainty on Tex reflects into a 
40-70% uncertainty on Miobe- 

One can estimate the outflow velocity from the terminal velocity of the CO spectra 
and the systemic velocity as Vaow = iVt — Vsys\- The dynamical time scale is then given by 
Ziobe/^ow, where the lobe length /lobe is defined as the maximum extent of the lobe measured 
from the 3 mm continuum peak. In all cases td lies in the range (1-3) xlO"^ yr. Note that 
this estimate of is to be taken as a lower limit because the lobe could extended beyond 
our field of view. 

The mass loss rate, Mfiow, was estimated from the ratio Miohe/td, and turns out to be 
on the order of 10~^ - 1O~^M0 yr~^ for all sources. Since molecular outflows appear to be 
momentum-driven (e.g., Cabrit & Bertout 1992) the momentum rate -Fco=M|obe^ow/^iobe 
may be taken as an indicator of the strength of the outflow and hence of the mass and 
luminosity of the YSO powering it. The G 16.59—0.05 outflow has F^o — 10^^ Mq km s~^ 
yr~^, whereas for G 28.87+0.07 and G 23.01—0.41 Fco are an order of magnitude higher. 
Values like these are typical of YSOs with luminosities of ~ 10'^ — 10^ Lq (see Figure 5 of 
Richer et al. 2000 and Figure 4 of Beuther et al. 2002a), confirming that we are indeed 
dealing with early B (proto) stars. 

These findings are not affected if one takes into account the unknown inclination angle 
i, here defined as the angle between the outflow axis and the line of sight. In fact, td, 
Mflow, and Fco, are proportional respectively to coti, tani, and sini/ cos^i. Thus, assuming 
as extreme values i = 20° and i = 70°, the corrections for the three quantities above are 
respectively 2.7-0.36, 0.36-2.7, and 0.39-8. In conclusion, for all reasonable inclination 
angles, the parameters listed in Table 7 might need a correction by less than an order of 
magnitude, thus leaving our conclusions unaffected. 

3.5. Temperature, Column Density and CH3CN Abundance of Cores 

3.5.1. Rotational Temperature and Column Density of the CH^CN Lines 

Symmetric top molecules like CH3CN are ideal probes to measure the gas kinetic tem- 
perature, because different K components belonging to the same J + 1 J transition are 
mainly excited by collisions, and spread over a few 10 MHz and can be observed simulta- 
neously in the same bandwidth (hence, suffering less calibration errors). In addition, for 
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densities typical of HMCs (> 10^~^ cm~^) the rotational temperature, Trot, obtained from K 
line ratios turns out to be very close to the kinetic temperature of the H2 gas. One may thus 
obtain reliable estimates of the temperature and column density of the region where CH3CN 
is detected through the rotation diagram method (Hollis 1982; Churchwell, Walmsley & 
Wood 1992). 

Since optical depth effects may affect the reliability of rotation diagrams, we have at- 
tempted to fit the K components of CH3CN with a simple LTE fit taking into account also 
the line opacity, as performed in e.g., Olimi et al. (1993). The results indicate that the 
optical depth is less than ~ 2, sufficiently low to allow usage of the rotation diagram method 
in our cases. 

To derive mean values of Trot and iVcHsCN over the cores using rotation diagrams, we 
analyzed the mean spectra in Figure [7] obtained by averaging the emission over the regions 
inside 5ct contour level of the K = 0+1 emission map. Assuming that all i^-components trace 
the same gas, we performed multiple- Gaussian fitting of the lines, forcing their widths to be 
identical and their separations in frequency to be equal to the laboratory values. Hence, the 
free parameters of the fit are the line intensities, the full width at half maximum (FWHM) 
of the lines, and the LSR velocity of one K component - arbitrarily chosen. Tables [8] and 
[9] summarize the fitting results, where Vsys is the best-fit LSR velocity. The line intensities 
were used to make rotation diagrams (see e.g. Churchwell, Walmsley & Wood 1992) from 
which we obtained an estimate of Trot and iVcHsCN- Note that our analysis is restricted to use 
the i^' = to 3 lines because of the limited correlator bandwidth. In these calculations, we 
assumed optically thin emission and adopted the partition function by Araya et al. (2005). 
The rotation diagrams for all three sources are presented in Figure [H where A^jk, -Ejk, and 
gjK represent the column density, energy, and statistical weight of level {J,K), and k is the 
Boltzmann constant. In Table [HI we give the obtained values of Trot and A'chsCN for the 3 
HMCs. 

The 3 HMCs have comparable mean A^chsCN, spanning the range (3-7) x 10^^ cm~^ 
(Table [9]), close to the values obtained from interferometric measurements in similar objects 
(e.g. G 24.78+0.08; Beltran et al. 2005). The mean values of Trot in G 16.59-0.05 and 
G 23.01-0.41 (120 - 130 K) are likely to be higher than that of G 28.87+0.07 (93 K), ahhough 
our Trot estimates have rather large errors. In G 28.87+0.07 both the temperature estimated 
from NH3 (see Table E) and the Tsb of the CH3CN lines (~ 5 K for = + 1) are twice 
smaller than in the other two sources (~ 10 K), suggesting that CH3CN gas should be colder 
than in G 16.59-0.05 and G 23.01-0.41. 

In Table [9|, one may notice that the mean Trot are a factor of 2 higher than the NH3 
kinetic temperature (T^) derived from the VLA observations (CTC97). One possible expla- 
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nation for this discrepancy is that optical depth effects in the CH3CN hnes may cause an 
overestimate of Tj-ot- It is also worth noting that the low-energy NH3 transitions observed 
by CTC97 may be affected by more extended emission (and hence colder gas) than those 
obtained through the CH3CN lines, as suggested in Olmi et al. (1993). In fact, these authors 
demonstrate that the two "thermometers", CH3CN and NH3, are in reasonable agreement 
for T ^ 50K, but for higher temperature the Trot estimated from CH3CN tends to exceed 
that from NH3. However, one should keep in mind that the relationship of Olmi et al. (1993) 
is based on single-dish observations with 6-8 times larger beam sizes than our interferometric 
measurements, and might hence be affected by material extended over larger regions than 
those imaged by us. 

3.5.2. Fractional Abundance of CH^CN Molecules in the Cores 

We have attempted an estimate of the fractional abundance of CH3CN (XchsCn) from 
the ratio between the mean A^^chsCN (TablelH]) and (A^(H2)) calculated in §3.1.2[ TableOgives 
the resultant XchsCN- We stress that the 3 mm emission from the three HMCs is dominated 
by thermal dust emission ( §3.1.21) . Note that the CH3CN line and 3 mm continuum emission 
show fairly similar spatial distributions (see panels (a) and (d) in Figures [9]- [TT]) . It is also 
worth pointing out that the errors on XchsCN take into account only uncertainties on the flux 
calibration (§ §2.1.11 and l2. 1.21) . One can see that XchsCN lies in the range (0.7 — 2) x 10~^, 
marginally lower than the values quoted for the well known HMC in Orion (5 x 10^^; see 
van Dishoeck et al. 1993). 

3.6. Velocity Structure and Stability of the CH3CN Cores 

Our results demonstrate that the ammonia cores imaged by CTC97 are indeed compact 
(~ 0.1 pc), dense (~ 10^ cm~^), and hot (~ 100 K) molecular cores. This strengthens the 
idea that high-mass (proto)stars - whose presence is suggested by the presence of H2O and 
OH masers (Foster & Caswell 1989) - are embedded in such massive (~ 100 Mq) HMCs. 
In the case of G 23.01—0.41, we also find evidence of a large scale clump seen in two CO 
isotopomers enshrouding the HMC, which is consistent with the findings of other authors 
(see e.g. Fontani et al. 2002 and references therein) in similar objects. In this section we wish 
to perform a more detailed analysis of the HMC structure by focusing on the velocity field 
of the gas, with the purpose of studying the HMC stability and infer the possible presence 
of rotation. 
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To gain insight into the velocity field of the gas, we used the same approach as Beltran 
et al. (2004), i.e. we determined the velocity in each point of the region where CH3CN is 
detected, by fitting the CH3CN spectrum in such points with the method described in §3.5.11 
In this way we also obtained the CH3CN line width {Aviat) in each point of the map: this 
is the FWHM after deconvolution of the instrumental resolution (Table [2]). In this process, 
we have considered only those points for which the intensity of the K = 3 line was above 
5a. The results are shown in Table [10] and Figures [9HTTJ It should be also noted that the 
regions over which the CH3CN velocity could be determined are only barely resolved, being 
comparable to the synthesized beam width. 

G 16.59— 0.06 The velocity field map of G 16.59— 0.05 (Figure [H^) shows a gradient ap- 
proximately along a line with P.A.~ 45°, roughly perpendicular to the outfiow axis. The 
most blueshifted gas lies at the center, while the most redshifted is seen to the northeast and 
southwest, close to the border of the region. A map of the line width is shown in Figure [9)d 
and appears to have no correlation with the velocity map. 

It would be surprising if the velocity gradient in this core were related to the CO outfiow 
shown in Figure HI because the two are inconsistent both geometrically and kinematically. 
On the other hand, it is difficult to believe that the CH3CN velocity is tracing rotation, since 
in this case the most blue- and redshifted gas should be found at the opposite extremes of 
the gradient. Furthermore, the velocity gradient corresponds to 1.2 km s~^ over 2400 AU, 
which implies an equilibrium mass of 0.2 Mq, much less than the mass of the HMC 

We believe, however, that the presence of rotation about the outfiow axis cannot be 
ruled out on the basis of our observations. In fact, the outfiow appears to be close to pole- 
on, which would make the component along the line-of-sight of the rotation velocity very 
small and hence difficult to measure. 

For a better understanding of the origin of the velocity gradient, one should also ex- 
plain the positional displacement between the 3 mm and CH3CN peaks. We propose two 
hypotheses for follow-up studies: 

• The continuum and CH3CN emission might arise from two distinct YSOs, the former 
deeply embedded in a relatively cold, massive core, the latter associated with hot gas, 
possibly distributed in a rotating ring, similar to those seen around binary systems in 
low-mass stars (e.g. GG Tau; Guilloteau et al. 1999); the lack of gas and dust at the 
center of the ring would explain why the continuum emission from the CH3CN peak is 
faint. 

• If the newly formed stars lie close to the surface of the core, this might cause a peak 



- 20 - 



in the CH3CN abundance close the the stars themselves and far from the core center, 
where the dust column density - and hence the mm continuum emission - peaks. In 
this scenario, the 3 mm peak could represent a high-mass protostar, while the CH3CN 
peak would correspond to a relatively more evolved OB star. 

G 28.87+0.08 The velocity map of the G 28.87+0.07 core (Figure ^) displays a clear 
gradient from southeast to northwest plus a less prominent one from northeast to south- 
west. It is worth reminding that the 3a contour levels of the 2.6 mm continuum (Figure [T)d) 
and CH3CN K = + 1 emission (Figure [5ll) are slightly elongated in a direction roughly 
perpendicular to the outflows ( §3.3.2p . Such an elongation is consistent with a flattened 
circumstellar toroid. 

The interpretation of the isovelocity map is complicated by the fact that the number of 
possible outflows and their structure is unclear in this region (see §3.3.21) . However, there 
are a few points worth of consideration: 

• The velocity gradients seen in CII3CN are inconsistent with the outflow velocity - no 
matter whether the one- or two-outflow hypothesis is correct - because the redshifted 
CH3CN gas is located approximately to the E-SE of the core, whereas the redshifted 
CO emission is found to the W-SW. 

• If a single outflow is present, its axis (solid line in Figure [TOb) runs perpendicular to 
the most prominent CH3CN gradient. 

• If instead two outflows are present, neither of these would be orthogonal to either 
CII3CN velocity gradient. 

Assuming that the CH3CN emission is tracing a core rotating about the axis outlined 
by a poorly coUimated CO outflow, the single outflow hypothesis in §3.3.21 seems preferable 
to the double outflow scenario. Additional evidence in favor of this scenario is given by the 
line width map of Figure [TUb. which shows that significant line broadening, by ~2kms~^, 
occurs from the core center to the border, just along the outflow axis. 

The CH3CN velocity gradient corresponds to a change of ~ 1.1 km s"^ over ~ 1.1 x 10^ 
AU. If this is due to rotation, the dynamical mass (Mdyn = -ResKot/^ with V^ot rotational 
velocity; Table [TD]) required to balance centrifugal and gravitational forces is equal to 3 ± 2 
Mq (Table [TOj) . This is an order of magnitude smaller than mass of the core estimated from 
3 mm continuum flux (Mj^^^ ~ 55 Mq] Table [TOj) . over the region where the rotation velocity 
is measured, and two orders of magnitude less than the virial mass (Mvij. ~ Rcs{Avij^ty / G ~ 
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1100 Mq; Table [TOj) obtained for the same region. This fact strongly suggests that even 
though rotation may be present in the HMC, its contribution to gravitational equilibrium is 
negligible. 

What does support the core against gravitational collapse? Here we should recall that 
the Af int map of Figure [TOb does not show a well-defined increase of the velocity width 
towards the core center, as one sees, instead, in similar objects - e.g., 031.41+0.31 and 
G 24.78+0.08 (Beltran et al. 2005; see their Figures 9 and 23). Therefore, our data do not 
provide any evidence for infall motions towards the core center, so that infall is unlikely to 
contribute to line broadening. A more likely explanation for the observed line width is that 
of turbulent motions of the core gas, as turbulence is believed to play a dominant role on 
the evolution of the core and the star formation process in its interiors (e.g., McKee & Tan 
2003). 

G 23.01—0.41 The G 23.01—0.41 core is significantly elongated in the southeast-northwest 
direction (see Figure [6li; P.A.~ 160°), namely perpendicular to the outflow axis, as one can 
see in the NH3, CH3CN, and 3.3 mm continuum maps in Figure [6l In the same figure, one 
can appreciate that also the 0^*^0(1-0) line presents a tail of emission in the same direction, 
extending towards the position of the 3 mm continuum peak. All the results strongly suggest 
that one could be observing a flattened structure rotating about the outflow axis. To shed 
light on this issue, as already done for the other two sources, we plot in Figure [TTh the velocity 
field of the CH3CN core. Regardless of some scatter, the velocity appears to increase quite 
steadily along the major axis of the HMG, as expected in case of rotation about the outflow 
axis (denoted by the solid line in the figure). The velocity shift amounts to 1.1 km s~^ over 
2.6 X 10^ AU. 

Unlike the velocity, the line width increases from north-east to south-west, roughly 
along the outflow axis. This suggests that to some extent the CH3CN emitting gas might 
also participate in the expansion. This fact could explain why the most blue- and redshifted 
CH3CN emission is observed respectively to the north and south of the HMC, although 
the main velocity trend lies along the SE-NW direction. We conclude, that, albeit slightly 
affected by the outflow, the CH3CN line emission from the HMC is mostly tracing rotation. 

As already done for G 28.87+0.07, we have estimated the mass needed for rotational 
support of the core. This is M^yn — 8 + 3 Mq, two orders of magnitude smaller than the 
HMC mass 

^d^rt^ - 340 Mq {Table M)- Instead, the virial mass is 5 times larger than 
the core mass, demonstrating that as well as for G 28.87+0.07, also for G 23.01—0.41 the 
dominant contribution to core equilibrium is not coming from rotation but from turbulence. 
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3.7. Comparison with the Other CH3CN Cores and Future Considerations 

As discussed in §3.6[ we believe that we have found significant evidence of rotation in 
two out of three cores, while in the third object detection of rotation might be hindered by 
projection effects (rotation axis almost parallel to the line of sight). In any case, rotational 
motions are undoubtly less prominent in these two cores than in those (G24.78+0.08 and 
G31.41+0.31) studied by Beltran et al. (2004, 2005). Given the low number statistics it is 
impossible to draw any conclusion out of this result. However, one cannot rule out the pos- 
sibility that this is an observational effect. The distance to G 23.01—0.41 is 1.4 times larger 
than that of G24.78 and G31.41, while the angular resolution of our observations is at least 
2 times worse than that of the PdBI observations by Beltran et al. (2004, 2005). Therefore, 
not only our observations are less sensitive to velocity gradients, but this makes it also more 
difficult to disentangle the contribution by distinct toroid-outfiow systems overlapping in the 
same region - especially when multiple sources are present as in the case of G 28.87+0.07. 
In conclusion, we believe that better angular and spectral resolutions are needed to make 
any search for rotating disks/toroids successful. 

Notwithstanding these caveats, it is worth stressing that the mere existence of rotation in 
two of the HMCs studied by us, albeit insufficient to guarantee support against gravitational 
forces, is an important clue for the process of high-mass star formation. In addition, one 
cannot rule out the possibility that, on smaller scales than those imaged by us, conservation 
of angular momentum might speed up the rotation, and thus attain centrifugal equilibrium 
in a circumstellar disk. This condition will be attained for radii satisfying the relation 



whose derivation is given in Appendix |Al A minimum value of -Roq can be obtained for 
p = 2.5 (see e.g. Fontani et al. 2002), which gives Req/Res — 6 x 10~^ and 5 x 10"^ for 
G 28.87+0.07 and G 23.01—0.41, respectively. Regions like these need sub-arcsec angular 
resolution and a higher density tracer than CH3CN to be investigated. 



Using the OVRO, Nobeyama, and IRAM-PdB mm-interferometers, we carried out in- 
tensive search for rotating toroids towards the massive YSOs in G 16.59—0.05, G 28.87+0.07, 
and G 23.01—0.41 which exhibit no or faint free-free emission (CTC97). Our observations 
revealed that these objects are embedded in HMCs with masses of 95-380 Mq and tempera- 
tures of 93-130 K, making the cores typical site of high-mass (proto)star formation. All the 




(1) 



4. 



Conclusions 
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3 objects harbored in the HMCs are driving powerful {Fco^ 10~^ — 10~^ MQkm yr~^) 
CO outflows. However, the nature of the outflows in G 28.87+0.07 and G 23.01—0.41 is un- 
clear; the origin of high velocity wing emission may be attributed to either single or double 
outflow(s). Such ambiguity made the interpretation of velocity gradients, identified through 
CH3CN i^-ladder line analysis, existing in the innermost densest part of the G 28.87+0.07 
and G 23.01—0.41 HMCs fairly difficult. The velocity gradients are almost perpendicular 
to their molecular outflow axes, suggesting the presence of rotating, flattened structures. 
However, the corresponding dynamical masses are an order of magnitude smaller than the 
masses derived from 3 mm dust continuum emission, which indicates turbulent pressure as 
the dominant support of the HMCs. No conclusion could be reached for the third source, 
G 16.59—0.05, as the putative rotation axis appears to lie close to the line-of-sight, thus 
making the detection of the rotation velocity very difficult for projection effects. Further 
higher resolution imaging will allow us to establish the presence of rotation on a more solid 
ground. 
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A. Centrifugal Radius 

Here we calculate an expression of the radius at which centrifugal equilibrium can be 
attained inside a rotating, spherically symmetric core. The hypothesis is that the rotation 
velocity V^ot observed at radius Rcs is not sufficient to sustain the core against gravitational 
forces thus allowing contraction and speed up of the gas in the core until centrifugal equi- 
librium is attained. If this occurs at radius Rcq with velocity V^q, conservation of angular 
momentum per unit mass gives the following relation: 



VeqReq — VrotReS 



(Al) 
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while the equilibrium condition turns into 

y^^.^GM^ (A2) 

The mass M{R) contained inside radius R depends on the dependence of the gas volume 
density n on the radius. Assuming that n varies with R~^, one obtaines M{R) oc 
Therefore, we can write the following relation 

M(i?,q) = M(/?efr) (l^y ' (A3) 

where M(i?efr) = M^^^^ measured by us from the mm continuum emission. 

Substituting Veq from Eq. (lAip and M(i?eq) from Eq. (]A3P into Eq. (]A2p . and replacing 
Vcot with Mdyn from the definition, we find 

= (A4) 
one consequently obtains the expression of i?eq as a function of measurable quantities: 
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Table 1. Summary of Interferometric Continuum Emission Observations 



Source 


Telescope 


Frequency 
(GHz) 


Bandwidth'' 
(GHz) 


SyiilliC'sizc't 

^maj ^ ^min 

(arcsec) 


1 JJt'aiii 
P.A. 
(deg) 


Image 

RMS 
(mjy beam~l ) 


G 16.59-0.05 


OVRO 


92.0'" 


3.5 


2.31 X 1.75 


-24 


0.50 




IRAM 


115.3 


0.32 


6.31 X 4.10 


-167 


2.49 


G 28.87+0.07 


OVRO 


92.0'" 


3.5 


2.10 X 1.80 


-74 


0.51 




IRAM 


115.3 


0.32 


4.97 X 4.15 


+23 


1.59 


G 23.01-0.41 


NMA 


98.8<= 


1.0 


2.69 X 1.93 


-14 


0.99 




NMA 


110.0<= 


1.0 


2.40 X 1.69 


-12 


1.35 




IRAM 


115.3 


3.5 


5.50 X 4.17 


-164 


1.90 



''Effective bandwidth of line- free channels for obtaining continuum image. 

''Center frequency of both the side bands. 
'^Center frequency of each single side band. 
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Table 2. Summary of Interferometric Molecular Line Observations 



Source 


Line 




Tracking 


Telescope 


Velocity 


Synthesized Beam 


Image'' 




Frequency 


Resolution 


^maj ^ ^min 


P.A. 


RMS 








(MHz) 




(km s-l ) 


(arcsec) 


(deg) 


(mjy beam~l ) 


G 16.59-0.05 


HCO+ 


(1-0) 


89188.518 


OVRO 


1.68 


3.26 X 1.95 


-10 


51 




CHsCN 


(5-4) 


90584.244 


OVRO 


1.63 


2.46 X 1.96 


-27 


18 




N2H+ 


(1-0) 


93173.900 


OVRO 


0.80 


2.45 X 1.85 


-46 


27 




12CO 


(1-0) 


115271.204 


IRAM 


1.63 


6.31 X 4.10 


-167 


43 


G 28.87+0.07 


HCO+ 


(1-0) 


89188.518 


OVRO 


1.68 


3.51 X 2.33 


-9 


23 




CH3CN 


(5-4) 


90584.244 


OVRO 


1.63 


2.47 X 2.05 


-61 


18 




N2H+ 


(1-0) 


93173.900 


OVRO 


0.80 


2.33 X 1.58 


-71 


28 




12CO 


(1-0) 


115271.204 


IRAM 


1.63 


4.58 X 4.58 


-21 


24 


G 23.01-0.41 


C180 


(1-0) 


109782.156 


NMA 


5.44 


7.44 X 5.52<: 


-21 


25 




HNCO 5()5 


-4()4 


109905.753 


NMA 


5.44 


7.44 X 5.52<: 


-21 


27 




CH3CN 


(6-5) 


109997.353 


NMA 


0.34'' 


2.86 X 1.67 


-9 


43 




13CO 


(1-0) 


110201.352 


NMA 


5.44 


7.44 X 5.52<= 


-21 


32 




12CO 


(1-0) 


115271.204 


IRAM 


0.41 


5.50 X 4.17 


-165 


43 



''Typical Image noise level with the velocity resolution shown in the column 5. 
''Effective resolution after smoothing. 

"^Received with the same intermediate-frequency, leading to the identical beam size for the lines. 
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Table 3. Properties of 3 mm Continuum Emission 



Source 




V 


Peak Position 


^cff 




(kpc) 


(GHz) 


R.A.(J2000) 


Dec.(J2000) 


(pc) 


(mjy) 


G 16.59-0.05 


4.7 


92.0 


18:21:09.06 


-14:31:47.9 


0.033 


37.5+7.5 






115.3 


18:21:09.06 


-14:31:47.5 


0.055 


79.6+16 


G 28.87+0.07 


7.4 


92.0 


18:43:46.25 


-03:35:29.9 


0.029 


11.3+2.3 






115.3 


18:43:46.31 


-03:35:29.2 


0.102 


42.6+8 


G 28.87+0. 07B 




92.0 


18:43:46.19 


-03:35:12.9 


0.033 


1.7+0.3 


G28.87+0.07C 




92.0 


18:43:46.31 


-03:35:15.4 


- 0.01 


1.5+0.4 


G28.87+0.07D 




115.3 


18:43:47.09 


-03:35:32.0 


0.030<i 


0.80+0.2 


G 23.01-0.41 


10.7 


98.8 


18:34:40.29 


-09:00:38.1 


0.068 


27.8+5.6 






110.0 


18:34:40.29 


-09:00:38.6 


0.055 


31.6+6.3 






115.3 


18:34:40.30 


-09:00:38.2 


0.103 


33.9+6.8 



Note. — See Table [T] for the parameters of observations and i|3.1.1l for the adopted detection 
criteria. 

''Distance to the region, see § 2 of CTC97. 
'^Effective radius, see i|3.1.1l for the definition. 

"^Total flux density integrated over the region enclosed by the 5(t level contour. 

^' \JA./tt for the area encompassed by the 5cr contour is given because the area is smaller than the 
beam size (Table [TJ. 



Table 4. Flux Densities of the 3 HMCs at Infrared Wavelengths 





21.3 fj,m 


14.7/jm 


12.1 




8.3 fim 


8.0^4111 


5.8 fj,m 


4.5 fim 


3.6 


fim 


Source 


S 




S 




S 




S 




5 




S 




S 




S 






(Jy) 


Note 


(Jy) 


Note 


(Jy) 


Note 


(Jy) 


Note 


(Jy) 


Note 


(Jy) 


Note 


(Jy) 


Note 


(Jy) 


Note 


G 16.59-0.05 


7.7 


1,2 


1.8 


1,2 


4.1 


1,2,3 


2.4 


1,2 


6.8E-2 


2,4 


3.3E-2 


2,4 


1.4E-2 


2,4 




5 


G 28.87+0.07 


14.7 


1,2 


6.8 


1,2 


5.1 


1,2 


2.8 


1,2 


1.3 


2,4 


8.0E-1 


2,4 


4.0E-1 


2,4 


6.0E-2 


2,4 


G 23.01-0.41 


4.4 


1,2 


< 0.51 


1,6 


< 0.86 


1,6 


< 0.20 


1,6 


l.lE-1 


2,4 


l.lE-1 


2,4 


6.3E-2 


2,4 


l.OE-2 


2,4 



Note. — 1: MSX data, 2: integrated inside the 50% level contour of the data, 3: contaminated with the nearby source(s), 4: GLIMPSE data, 5: impossible to define a point source (see Figure 
[T] and 6: 3(7 upper limit. See i|3. 1.21 for the details. 
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Table 5. Core Mass Estimated from 3 mm Continuum Emission 



Source 


Mduet" 


(7V(H2)>'^ 


(Mq) 


(xlO^^ cm-2) 


G 16.59-0.05 


95±6 


3.0±0.9 


G 28.87+0.07 


100±7 


4.2±1.3 


G 23.01-0.41 


380±22 


3.6±1.1 



Note. — See i|3.1.2l for the details. 

''Mass of the core estimated from the 3 mm 
continuum flux density shown in Table |3] We 
assumed that T;just is equal to Trot of CH3CN 
emission (Table [9]l. 

''Mean H2 column density calculated over 
the area encompassing the emission inside the 
5(T contour. 



Table 6. Properties of the Large-Scale Clumps Surrounding the 3 HMCs 











Source 


Probe 


(pc) 


(A^q) 


G 16.59-0.05 


HCO+ (1-0) 


0.070 


440 - 1500 




N2H+ (1-0) 


0.15 


540 


G 28.87+0.07 


HCO+ (1-0) 


0.12 


84 - 360 


G 23.01-0.41 


i^CO (1-0) 


0.426 


1 X 10-* 




CISQ (1-0) 


0.174 


2 X 10"' 



Note. — For details, see 

''Effective radius, see i|3.1.1l for the definition. 

^ LTE-mass obtained by integrated the emission over the 
area enclosed by the 5cr level contours. 
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Table 7. Parameters of Molecular Outflows 



Source 


Lobe 




id"- 












(Mq) 


(yrs) 


{Mq yr-i) 


{Mq km s"-"- yr" 




G 16.59-0.05 


Blue 


21±11 


2.7x10-* 


(2±l)xl0-* 


{9±5)xl0-3 






Red 


7±4 


1.4x10* 


(5±3)xl0-* 


(9±5)xl0-3 




G 28.87+0.07 


Blue 


25±12 


2.6x10* 


(3±2)xl0-* 


(2±l)xl0-2 






Red 


90±40 


3.1x10* 


(2±l)xl0-3 


(5±3)xl0-2 




G 23.01-0.41 


Blue SW 


17±12 


2.5x10* 


(9±7)xl0-* 


(2.5±1.5)xl0- 


■2 




Blue NE 


14±10 


2.8x10* 


(6±3)xlO-* 


(1.9±1.2)xl0- 


■2 




Red SW 


140±70 


3.4x10* 


(4±3)xl0-3 


(1.3±0.7)xl0- 


-1 




Red NE 


4±3 


3.4x10* 


(1.2±0.7)xl0~* 


(3.5±3.0)xl0- 


3 



Note. — Outflow characteristic estimated from the *^CO (1-0) data, see tj3.4l All the errors 
are considered only for the Tox uncertainty. 

''Outflow lobe mass. 

''Dynamical time scale. 

'^Outflow mass loss rate. 

Outflow momentum flux. 



Table 8. Results of CH3CN Spectra Analysis 



Source 


Transition 


VLsr" 






/ T^bdv (K-km s-i) 




(km s-i) 


(km s~^) 


K = 


1 


2 


3 


G 16.59-0.05 


J =5-4 


59.9+0.2 


5.6+0.2 


20.5+1.2 


14.1+1.2 


10.5+1.1 


13.7+1.1 


G 28.87+0.07 


J =5-4 


103.5+0.3 


9.1+0.4 


12.9+1.4 


11.5+1.5 


7.8+0.9 


7.7+0.9 


G 23.01-0.41 


J =6-5 


77.4+0.1 


8.3+0.2 


16.2+0.8 


16.7+0.8 


15.0+0.6 


13.8+0.6 



Note. — For details, see i|3.5.1l 

''Centroid velocity from the multiple-Gaussian profile fitting. We adopted the Vlsr in column 3 as the Vs 
of the HMCs. 

^ Line width in FWHM. 



Table 9. Summary of Methyl Cyanide RD Analysis 







Temperature'' (K) 






Source 










(pc) 


T,ot(CH3CN) Tk{NH3) 


{cm-2) 


(xlO-iO) 



G 16.59-0.05 0.029 ISOtas 54 (6.9lig)xl0" 20±7 

G 28.87+0.07 0.042 93^11 37 (3.ltJ|)xlO" 7±2 

G 23.01-0.41 0.065 ISltlg 58 (4.6tog)xlO" 13±4 



Note. — For details, see ^ i|3.5.1l and l3.5.2l 

''Effective radius for the CH3CN K = + 1 emitting region, see i|3.1.1l for the definition. 
''Rotational temperature for CH3CN ( ^3.5.11 1 and kinetic temperature for NH3 (CTC97). 
''Total column density of CH3CN. 

'^Fractional abundance of CH3CN obtained from comparisons with the mean (A'^(H2)) in Table[5] 



Table 10. Properties of the Central Regions of the 3 HMCs 



Source 








dust 




(Mq) (km s-i) 




(km s-l) (Mq) 





G 16.59-0.05 95±6 5.4±0.3 180±15 

G 28.87+0.07 100+7 9.1+0.4 780+78 ~ 0.5 3+2 55+3 

G 23.01-0.41 380+22 8.6+0.2 1100+40 ~ 0.6 8+3 340+20 



Note. — For details, see ^3^ 
''From Table [5] 

''Intrinsic velocity width (Aiiint) in FWHM after deconvolving the instrumental ve- 
locity resolution (Table [2]l where Anint is measured over the region enclosed by the 
50% level contour of the emission. Here, contributions of thermal gas motions are 0.38, 
0.32, and 0.37 km s"' for G 16.59-0.05, G 28.87+0.07, and G 23.01-0.41, respectively, 
when we assume that kinomatical temperature (T]j) of the gas is equal to the mean Trot 
in Table [9] 

''Virial mass over the 50% area of CH3CN K = + 1 emission. We adopted ij^ff in 
Table |9] for calculating M^ir- 

Rotation velocity assuming that the velocity gradient is produced by rotation of the 
core. 

''Dynamical mass. 

'Core mass calculated from the 3 mm continuum flux density over the region enclosed 
by the 5a level contour of the CH3CN K = 3 emission. 
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(a) G16.59-0.05 



3 1 ^m+i.3mni T1^zm+5.Gmm S.Opru l-Gfun 




1^21 0H*+ 

(b) G28,87+0,07 



3 r >;m+3.3mni ?1>Tn+2.amm a.Ptim 3.6pin 




1 1 



Fig. 1. — Continuum emission maps towards the 3 HMCs at mm- and infrared wavelengths. The 
left two panels for each source shows an overlay of mm-continuum emission (white contours) on 
the 21 /im one from the MSX satellite data base; the single white dashed contour indicates the 
50% level of the 21 /im emission. The ellipses in the inserted panels show the HPBWs of the 
synthesized beam (Tabled]). The right two panels show 8.0 and 3.6 /im images from the Spitzer 
GLIMPSE survey; the stars mark the peak positions of the mm continuum source, and the thin 
contour indicate the 50% level contour of the peak intensity of the main cores of our interests. 
Contours for the mm-continuum maps are 2a intervals, starting from the 3a level (see Table [T] for 
the RMS noise levels). In the panel (b), the associated labels indicate names of the mm sources 
identified by us (see i j3.1.1l and Table [3]) . 
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Fig. 2. — Continuum spectra of the 3 HMCs. The flux densities at 3 mm and infrared bands 
are summarized in Tables [3] and HJ respectively. 
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Fig. 3. — Interferometric molecular line spectra towards the peak positions of each emission in 
Tsb scale. The black vertical dashed lines indicate the Vsys obtained from the CH3CN line analysis 
(TableED: Vlsr = 59.8 km s'^ for G 16.59-0.05, 103.3 km s'^ for G28.87+0.07, and 77.3 km s"^ 
for G 23.01—0.41. The vertical blue and red dashed lines throughout the panels show the boundary 
velocities (14>; see ^3.3.1l for the definitions) between the core gas and the blue- or redshifted outflow 
lobes, respectively. The blue, green, and red bars under the spectra indicate velocity ranges to 
obtain integrated intensity maps shown in Figures U] - [6j The CH3CN emission seen in l^sR ^ 60 
km s~^ at G 23.01—0.41 is the other ii'-components (see Figure [7]). 
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Fig. 4. — Maps of molecular lines and continuum emission for G 16.59—0.05. (a) Overlays of the 
blue- and redshifted wing emission maps of ^^CO (1-0) on the 3 mm continuum emission (grey 
scale), total integrated intensity maps of (b) HCO"^ (1~0) (c) N2H^ (1^0) emission (grey scale 
plus contour), and (d) overlay of total maps of NH3 (3,3) (green contours; CTC97) on CH3CN K = 
and 1 line emission (grey scale). In the panels (b)-(d), the thin blue- and red contours indicate 
the 5a levels of the blue- and redshifted ^^CO (1-0) lobes shown in (a). All the contours, except for 
the ^^CO, start from the 3a level with a 3a step. The dashed contours are correspond to —3a and 
-6a levels for HCO+ in panel (b), N2H+ in (c), and CH3CN in (d). Contours for the ^^CO outflow 
maps are drawn with 5a, 10a, 15a, 20a, and 40ct levels for the clarity of the plots. The RMS noise 
levels of the images are 33, 19, 0.50, 0.43, 0.10, 4.2, and 0.26 mJy beam"! for ^^CO blueshifted, 
^■^CO redshifted, 3mm continuum, HCO"*", N2H+, NH3, and CH3CN maps, respectively. The solid 
lines in (a) indicate the identified outfiow axes ( ^3.3.ip . The central star in (b)-(d) marks the 
peak position of the 3 mm continuum emission. The blue and redshifted ^^CO (1-0) emission is 
integrated over the velocity ranges of 47.4 < VLsR/kni s^^ < 54.0 and 63.4 < VLSR/km s^^ < 78.1, 
respectively. The HCO"*" (1-0) line is integrated over 53.3 < VcsR/km s~^ < 66.7. These velocity 
ranges are shown by the horizontal color bars under the spectra shown in Figure [3^. For the N2H"'" 
(1-0) emission, all the hyperfine emission detected in 49.4 < V^gR/km s~^ < 70.6 is integrated. 
The frequency range to obtain the CH3CN K = + 1 map is shown by a horizontal bar under 
the corresponding spectrum in Figure [71 The ellipse at the lower corners in each panel shows the 
synthesized beam size (Tables [U and [2]) . 
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Fig. 5. — Maps of molecular lines and continuum emission for G 28.87+0.07. All the contours and 
symbols are the same as those in Figured! The RMS noise levels of the images are 47, 24, 0.51, 
0.3, 0.13, 5.0, and 0.18 for ^^CO blueshifted, ^^CO redshifted, 3mm continuum, HCO+, N2H+, 
NH3 (CTC97), and CH3CN maps, respectively. The dashed contours correspond to the —3a level 
for HCO+ in panel (b), N2H+ in (c), and CH3CN in (d). The blue and redshifted ^^CO (1-0) 
emission, HCO"*", and N2H"'" are integrated over the velocity ranges of 83.3 < VLSR/km s^^ < 96.3, 
111.8 < lisR/km s-i < 123.9, 102.3 < VLsn/km s'^ < 114.0, and 95.0 < VLSn/km s~i < 111.1, 
respectively (see also Figure [3)3). The star with the associated labels indicate the peak positions 
and names of the mm sources identified by us (see ^3.1.11 and Table [3]). All the other symbols are 
the same as in Figure HI 
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Fig. 6. — Maps of molecular lines and continuum emission for G 23.01—0.41: (a) overlay of blue- 
and redshifted wing emission maps of ^^CO (1-0) on the 3mm continuum emission (grey scale), 
(b) overlay of blue and redshifted wing emission maps of the ^^CO (1-0), total integrated intensity 
maps of NH3 (2,2) (green contour; CTC97) on the integrated intensity map of the bulk ^^CO (1-0) 
emission (grey scale plus contour; 69.0 < ^LSR/km < 85.3), (c) total integrated intensity map 
of C^^O (1-0) (grey scale; 60.0 < ^LSR/km s~^ < 87.2), and (d) overlays of the blue and redshifted 
HNCO 5o5 — 4o4 emission on the total map of CH3CN (6-5) K = + 1 emission (grey scale). 
The dashed contours correspond to the —3a level for ^'^CO in panel (b), C^^O in (c), and HNCO 
in (d). The blue- and redshifted wing emission are integrated over the velocity ranges of 49.9 < 
FLSR/km s-i < 71.5 and 86.0 < VLsn/km s'^ < 101.6 for ^^CO, 47.2 < VLsn/km s"^ < 69.0 and 
85.3 < ^LSR/km s-i < 96.2 for ^^CO, 48.1 < VtsR/km s"^ < 69.9 and 86.2 < VLSR/km 3"^ < 
91.6 for HNCO, respectively (see also Figure [3h). All the contours and symbols are the same as 
in Figure H The RMS noise levels of the images are 90, 26, 1.4, 0.70, 0.79, 5.1, 21, 22 and 0.13 
mJy beam-i for ^^CO blueshifted, ^^CO redshifted, 3mm continuum, ^^CO, C^^O, NH3, HNCO 
blueshifted, HNCO redshifted, and CH3CN maps, respectively. All the other symbols are the same 
as in Figure m 
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Fig. 7. — Mean CH3CN spectra integrated over the 5cr level areas of the 3 HMCs in flux density 
(S^) scale: J = 5-4 for G 16.59-0.05 and G 28.87+0.07, and J = 6-5 for G 23.01-0.41. The 
vertical bars above the spectra indicate the rest-frequency of the K rotational ladder emission. 
The horizontal bars under K = and 1 emission indicate frequency ranges to obtain the integrated 
intensity maps in the panels (d) of Figures U] - [6l 



-42- 



26 



1 \ r 
(a) 



^ ^ \ ' 
G1 6.59-0.05 



25 



1" 



24 



_T,„, = 130 K 
N(CH3CN) = 6.9E+14 cm" 

_J \ \ \ \ 



J L 



26 



25 



24 



_ 1 1 1 


1 1 1 _ 




G28.87+0.07 _ 




CH3CN 5k-4k 


_ = 93 K 




N(CH3CN) = 3.1E+14 


-2 

cm 


-III 


III- 



40 



EjK/k (K) 



26 



80 



40 80 
EjK/k (K) 



25 



24 



_ 1 1 1 


1 1 1 _ 




G23.01-0.41 _ 




CH3CN 6k-5k 




i — 






_Tro, = 121 K 




N(CH3CN) = 4.6E+14 


-2 

cm 


-III 


III- 



40 



EjK/k (K) 



80 



Fig. 8. — Rotational diagrams for CH3CN Jk — ( J — 1)k transitions obtained from the 
averaged spectra over the 5a areas of the = + 1 emission (Figure [7]), see §3.5.11 The 
natural logarithm of the column density in state J per substate [In^Nj^/ g^x)] (cm-^) is 
plotted to the upper energy level of Ej/k. The dashed straight line is a least-square fit to 
the data, yielding rotational temperature (Trot) and column density, iV(CH3CN), shown in 
each panel. All the derived parameters are summarized in Table M 
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Fig. 9. — (a) Isovelocity and (b) intrinsic velocity width (A^int) maps obtained from the multiple 
Gaussian line-profile fitting to the CH3CN emission towards the central portion of the G 16.59—0.05 
HMC. These maps are presented inside the 5a level contour (inner solid contour) of the K = 3 
total maps; the outer solid contour indicates the 5a level of the K = + 1 map. The dashed 
thin contours show the positive contour for the K = 0+1 map in Figure [H The ellipses at the 
bottom-right corners indicate synthesized beam size (Tabled]). The straight lines and central stars 
are the same as in Figure H Systemic velocity (Kys) of the HMC is Vlsk= 59.8 km s"^ (Table ED, 
and thermal line width at = Trot = 130 K is 0.38 km s~^. 
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Fig. 10.— Same as Figure El but for G 28.87+0.07. The V^ys is Vlsr= 103.5 km s'^ (Table ED, 
and thermal line width at 93 K is 0.32 km s~^. 
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Fig. 11. — Same as Figure[9l but for G 23.01—0.41. The transition employed here is J =6~5. The 
Vsys is Vlsr= 77.4 km (Table [8]), and thermal line width at 121 K is 0.37 km s~^. 



